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Abstract 

The Richardson exact solution for the reduced BCS Hamiltonian is applied to examine 
how sensitive are the pairing characteristics (condensation energy, spectroscopic gap, parity 
gap) to a specific configuration of single-electron energy levels in nanosize metallic grains. 
Using single-electron energy spectra in parallelepiped-shaped potential boxes with various 
volumes and aspect ratios as a model of energy levels in grains, we show that this sensitivity 
is extremely high. Just due to such an extreme sensitivity, the detailed shape of grains 
cannot be detected through the pairing characteristics, averaged over an ensemble of grains, 
even in the case of relatively small size dispersion within this ensemble. We analyse the 
effect of the pairing interaction on the excited-level spacings in superconducting grains and 
comment on the influence of shape- dependent fluctuations in single-electron energy spectra 
on the possibility to reveal this effect through tunnelling measurements. 



I. INTRODUCTION 

The experiments of Ralph, Black and Tinkham [0-0 , revealing, in particular, the pres- 
ence of a parity-dependent spectroscopic gap in tunnelling spectra of nanometre-scale Al 
grains, inspired intensive theoretical studies of pairing correlations in those grains. Recently, 
the exact solution to the reduced BCS Hamiltonian, which was developed by Richardson 
in the context of nuclear physics as long ago as in 1963 [|J, has been reintroduced to the 
condensed matter community and used to describe superconductivity in nanosize metallic 
grains by Brown and von Delft || . It is noteworthy that the Richardson solution is applica- 
ble at arbitrary distributions of single-electron energy levels, which form, together with the 
interaction strength, a complete set of input parameters for the reduced BCS Hamiltonian. 
This allows one to go beyond the model with equally spaced energy levels that is often used 
for nanosize superconducting grains. Thus, the Richardson solution has been applied || 
to study the effect of level statistics on superconductivity in nanosize grains with random 
levels, assumed to follow the Gaussian orthogonal ensemble distribution, and to check the 
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quality of the previous treatment of this effect 0. Keeping in mind that an appropriate 
choice of statistics to describe the energy-level distribution in concrete nanostructures might 
be not self-evident, it seems interesting to analyse in more detail how sensitive are the pair- 
ing characteristics to a specific configuration of "particle-in-a-box-like" || single-electron 
states in grain. Based on the Richardson solution [|J, in the present communication we 
tackle this problem by modelling single-electron energy spectra in grains with those spectra 
in parallelepiped-shaped potential boxes of various shape and size. The influence of the 
grain shape on the pairing characteristics (condensation energy, spectroscopic gap, parity 
gap) is studied both for single grains and for their ensembles. We also consider the effect 
of the pairing interaction on the excited-level spacings in superconducting grains and touch 
upon the influence of shape- dependent fluctuations of single-electron energy levels on the 
possibility to detect this effect through tunnelling measurements. 



II. THE MODEL 

For a nanosize grain, the BCS pairing Hamiltonian can be written as (see, e. g., |8|||) 

H = E s i a )^j<r - Ad E a \Al a n a ih (!) 

where the operator aj CT (dj a ) creates (annihilates) an electron in the j-th time-reversed 
single-electron state with the spin o and the energy Sj. The second term in Eq. (HD is the 
interaction Hamiltonian. The sum in this term is over the set I of Ji states inside the energy 
interval [e-p — Tuod^f + hui£>), which will be referred to as the interaction band, with the 



Debye energy hup = 34 meV for Al [1(J and £p, the Fermi energy. The interaction strength 
is a product of the mean energy-level spacing within the interaction band, d = 2huu/Ji, 
and the dimensionless parameter A, taken to be 0.22 (close to that for Al |TT|). 



Since the electrons occupying levels outside the interaction band are straightforwardly 
described by the first term in the r.h.s. of Eq. (p]), we will focus only on the electrons, 
which reside in the interaction band. Let the interaction band be populated by N = 2n + b 
electrons, of which b electrons are on singly occupied levels (the set of these levels, blocked 
to the pair scattering, is denoted as B). The electrons on singly occupied levels do not 
participate in the pair scattering and their contribution to the energy of the iV-electron 
system under consideration is simply J2jeB e j- The remaining 2n electrons form n pairs, 
which are distributed among the set U = I \ B of Jj — b unblocked levels. Richardson Q 
showed that the energy of these pairs is given by the sum of n parameters E\, . . . , E n , which 
are non-degenerate roots (E^ ^ E v for all // ^ v) of the following set of n coupled equations: 

1 ^ ~\ d ^ ^ -\ 2c/ q (2^ 

A j eU 2ej — E v M=1 E^ — E u 

v — 1, . . . , n . 

Thus the total energy of electrons, which populate the interaction band, is 

n 

E N ,b = E £ J + E E», N = 2n + b. (3) 
jeB v=l 
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Given the energies £?jv,6, the pairing characteristics like the condensation energy, the 
spectroscopic gap, and the parity gap can be calculated. Our first aim is to analyse how 
sensitive are these characteristics to a specific configuration of single-electron energy levels, 
Ej, in the interaction band of a grain. As an example, we will consider the single-electron 
energy spectra in a parallelepiped-shaped hard-wall box with sizes l X i ly cincl lz- 

„ _ 7T 2 ft 2 (nl(j) nl(j) nl(j)\ 

J ~ 2m { 11 + 11 + 11 ' [ ] 



n x {j),n y {j),n z (j) = 1,2,3 



where m is the effective mass of an electron (m = 1.4m e in Al HlQf). The correspondence 
between j and n x (j), n y (j), n z (j) is determined by the requirement Ej < £j+i, which we 
impose. When labelling the single-electron energy levels, we will take j = for the upper- 
most energy level populated by electrons in the ground state. A grain with sizes l X i ly and l z 
contains N tot = n e l x l y l z electrons, where the electron density is taken to be n e = 181 nm -3 
for Al [|Hj. A major part of these electrons are on energy levels below the interaction band, 
which are assumed to be completely filled. Assuming in addition that all the energy levels 
with 6j > Ep + hu>D are empty, the number of electrons in the interaction band, N, is fully 
determined by the volume and the shape of the grain, the parity of N being the same as 
that of iVtot- 



III. SHAPE-DEPENDENT PAIRING CHARACTERISTICS OF A NANOSIZE 

GRAIN 

First, we will consider the condensation energy, which can be expressed as (see, e. g., 
Refs. 

E^ b (X) = E% b (0)-El b (X)-nXd, (5) 

where E^ b is the ground-state energy of the A-electron system in the interaction band 
(6 = for even N = 2n and 6=1 for odd N = 2n + 1). In Fig. 1 we show the calculated 
condensation energy for nearly cubic grains with a fixed number of electrons, A tot = 4000 
(the corresponding volume of grains is approximately 22 nm 3 ), when varying the aspect 
ratio, l x : l y : l z — 1 : (1 + rf) : (1 + 2rf). As seen from Fig. 1, the condensation energy 
exhibits well-pronounced maximums at some values of the parameter 77. An inspection of 
the corresponding single-electron energy spectra, Ej, shows that the maximums appear when 
a degeneracy or quasi- degeneracy occurs between the lowest empty energy levels and the 
uppermost occupied energy levels. The degree of degeneracy is especially strong in a perfect 
cube resulting in a huge (about 104 meV high) peak of the condensation energy at rj = 0. 



However, this peak is seen to be extremely narrow implying that a proposal |T3[ for an 
experimental observation of degeneracy effects typical for perfectly symmetric grains can 
be hardly realised. At the same time, even for less symmetric grains, the variations of the 
condensation energy as a function of the grain shape are still strong (about one order of 
magnitude) reflecting a high sensitivity of the pairing characteristics to the detailed structure 
of the single-electron energy spectrum in the vicinity of the Fermi level. As illustrated by 
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Fig. 2, this sensitivity is also manifested when considering the condensation energy for grains 
with a fixed aspect ratio as a function of their volume or, equivalently, of N tot . 
Figure 2 shows also the lowest excitation energies, 

{ . } _ fE$ tl -E% tl for N = 2n + 1, 

N ~\E%-E- Nfi foriV = 2n, (6) 

z = l,2,3,..., 

relevant to the so-called "primary" peaks fT2j in tunnelling spectra of nanosize supercon- 
ducting grains. In grains with an odd number of electrons, the energies i%\ correspond 
to excited states with a single unpaired electron on the i-th level (i > 0). For a grain 
with an even number of electrons, the energies E$ 2 correspond to excited states, where two 
electrons are unpaired, one of them having the energy e an d another occupying the i-th 
single-electron energy level. Since these excited states involve one broken Cooper pair, the 
differences between their energies and the energy of the fully paired ground state are at 
least as large as the energy cost for pair breaking. The quantity d$ defined for an even 
N coincides with the spectroscopic gap |14|]. Figure 2 allows one to see a close similarity 
between the behaviour of the condensation energy as a function of N tot and that of the 
contribution to the spectroscopic gap due to the pairing interaction. Another observation, 
implied from Fig. 2, is an increase of the density of excited energy levels in a superconduct- 
ing grain as compared to that in a normal grain (with A = 0). The physical reason for this 
increase is the blocking of energy levels by unpaired electrons that reduces the interaction 
energy of the remaining electron pairs. With increasing i, the effect of the i-th level blocking 
on the interaction energy relaxes. Therefore, a downward energy shift due to the pairing 
interaction is larger for higher excited states than that for lower ones. This means that the 
energy spacing between excited states decreases when the interaction is switched on. Such 
a decrease was recently analysed within a perturbative approach for highly excited states 
(si — 6Q^>d) in the case of equidistant single-electron energy levels [12]. However, it follows 
from the above discussion as well as from Fig. 2 that an increase of the energy-level density 
due to the pairing interaction is more pronounced just for the lowest excited states and for 
small inter-level spacings £j — Eq, when a perturbative treatment is hardly applicable. 
In Fig. 3, we plot the calculated parity gap 0, 

= #2^+1,1 — - (^2^,0 + ^2n+2,o) > (?) 

i. e. the difference between the ground-state energy of a grain, which contains an odd total 
number of electrons, 2n+l of them being in the interaction band, and the mean ground-state 
energy of the grains obtained by removing or adding one electron. Since Eq. (0) involves the 
ground-state energies of three systems with different number of electrons, the behaviour of 
as a function of A^ot is somewhat richer as compared to that of the condensation energy, 
E 2n+ i i, which is also shown in Fig. 3. 



IV. PAIRING CHARACTERISTICS FOR ENSEMBLES OF NANOSIZE GRAINS 

Having demonstrated an extreme sensitivity of pairing characteristics to the specific 
geometry of a single grain, we will now analyse these characteristics for ensembles of nano- 
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size grains. Let us consider a set of parallelepiped-shaped grains whose sizes are given by 
expressions 

l a = (L) + r a Al, a = x,y,z (8) 

with random values r a e (—1, 1). Since we will deal with values Al < 0.25 nm, which are 
expected to be much smaller than the size dispersion for ensembles of realistic grains, all 
values of r a within the interval (—1, 1) can be taken equally probable. For several ensembles 
of 500 to 700 grains with (l x ) : (l y ) : (l z ) = 1:2:2 and Al = 0.2 nm, the ensemble average 



\2 



1/2 

are 



for the condensation energy, (E^ ), and the variance SE^ = ((_E^ ) 2 ) — (E^ ) 
shown in Fig. 4(a) as a function of (l x ). Analogous plots for the spectroscopic gap and its 
variance are given in Fig. 4(b). For comparison, we display the condensation energy and the 
spectroscopic gap in grains, where single-electron energy levels in the interaction band are 
equally spaced by (d), the average energy distance between adjacent levels for corresponding 
ensembles. For grains with equidistant levels Ej, the condensation energy is appreciably 
smaller than (E^ ). As concerns the spectroscopic gap, the difference between the results 
for grains with equidistant levels and those for ensembles of grains is less pronounced. There 
are two (positive) contributions to the spectroscopic gap: one is due to the energy spacing 
E\ — Sq, another is determined by the coupling energy of an electron pair. Since the coupling 
energy increases with decreasing E\ — Eq, fluctuations of these two contributions within an 
ensemble of grains partially cancel each other. 

It is worth noting that for the ensembles of grains under consideration, the variances 
SE^ and 5d$ as well as the differences between the ensemble averages (E^ ) ((djy )) and 

the corresponding values E^ (d$) in grains with equidistant levels are several times larger 
as compared to the results of Ref. ||, where single-electron energy levels in grains were 



assumed to follow the Gaussian orthogonal ensemble (GOE) distribution [15]. The reason 
for this discrepancy is that for the single-electron energy spectra given by Eq. (J|), energy- 
level spacings obey the Poisson distribution. In particular, small nearest-neighbour spacings 
(sj — £j_i <C d), which are unfeasible in GOE, have the maximum probability density in 
the ensembles analysed here. It is interesting to address the question of whether or not the 
"average grain shape", which is given in our example by the aspect ratio (l x ) : (l y ) : (l z ), can 
be revealed through the pairing characteristics of an ensemble of grains. Since the variances 
5Eft and Sd^ are relatively large, the answer appears to be negative. Thus, the relative 
differences between the averages, calculated for three ensembles of grains with approximately 
one and the same volume {(l x )(ly)(lz) — 172.5 nm 3 , Al = 0.2 nm) but substantially different 
aspect ratios ((l x ) : (l y ) : (l z ) = 1 : 1 : 1, 1 : 2 : 2, and 1:2: 16), do not exceed 2.5% for 
(E% fi ) and 3% for 

Let us analyse the aforementioned effect of the pairing interaction on the spacings be- 
tween excited states with energies E$ b (b=l for odd N and 6 = 2 for even N). Counting 

the energy from the lowest excited level, E^\, we can write down the density of higher 
energy levels as 

D(E)=Y,S(E-[E%-E$ b ]). (9) 

i>2 

Since the energy spectra are strongly affected even by weak fluctuations of the grain shape, 
modifications of the "detailed" level density D(E) due to the pairing interaction might be 
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hidden by seemingly irregular distribution of energy levels, especially when studying a small 
grain with only a few energy levels in the interaction band. As a more convenient quantity, 
we introduce the excited- level density averaged over a relatively large energy interval A: 

([B/A]+1)A 

D A (E) = - / D(E)dE, (10) 

[E/A]A 

where [x] denotes the integer part of x. To reveal more clearly the influence of the pairing 
interaction on the excited-level density, we examine the function 

D A (E)-D £ (0) 

PaAE) ^ , (11) 

which describes relative deviations of D&(E) from the level density averaged over an energy 
interval £ 3> A. In Fig. 5, we show the ensemble averages, (pa,s), calculated for two en- 
sembles of 1000 nanosize grains ((l x ){l y ){lz) = 23.3 nm 3 ) with an even [panel(a)] and odd 
[panel(b)] number of electrons. In spite of fluctuations of (pa,£), which are still present even 
for relatively large ensembles under consideration, a pronounced trend to an increase of the 
excited-level density when lowering E is seen from plots for superconducting (A ^ 0) grains 
at any parity of the number of electrons. In this respect, the behaviour of (pa,e) for super- 
conducting grains is quite distinctive from that for corresponding normal (A = 0) grains. 
The inset in Fig. 5(b) demonstrates that for largish grains (with (l x )(ly)(lz) = 256 nm 3 ) 
the aforementioned distinction can be revealed even when averaging the function p^,£ over 
a few-grain set. The obtained results imply that the excited- level densening due to the 
pairing interaction might be detected through experimental tunnelling spectra accumulated 
over sets of grains. 
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Figure captions 



Fig. 1. Condensation energy E^ Q , calculated according to Eq. (||) for parallelepiped-shaped 
grains with the aspect ratio l x : l y : l z = 1 : (1+??) : (1 + 27/), as a function of the parameter rj. 
The volume of grains is fixed to keep a constant number of electrons in a grain, iV to t = 4000. 
The number of electrons in the interaction band, N, varies from 20 to 42 for the range of rj 
under consideration. 

Fig. 2. Condensation energy E^ , given by Eq. (|5|), and the lowest excitation energies 

d$ (i = 1 to 5), given by Eq. (|6]), for parallelepiped-shaped grains with the aspect ratio 
lx '■ ly '■ lz = 1 : 1-3 : 1.7 as a function of the number of electrons iVt ot . Lines are guides for eye. 
Only grains with even iV tot are considered. The number of electrons in the interaction band, 
N, varies from 108 to 146 depending on N tot . For comparison, the five lowest excitation 
energies for the same grains in their normal state (at A = 0) are also shown. 

Fig. 3. Parity gap A^, given by Eq. (|7]), and condensation energy -E , ^ +1 1( given by Eq. (0), 
for parallelepiped-shaped grains with the aspect ratio l x : l y : l z = 1 : 1.9 : 2.1 as a function 
of the (odd) number of electrons iVtot- Lines are guides for eye. The number of pairs in the 
interaction band, n, varies from 83 to 110 for grains under consideration. 

Fig. 4. Ensemble averages and variances of the condensation energy [panel (a)] and of the 
spectroscopic gap [panel (b)] as a function of the average grain size (l x ) for ensembles of 500 
to 700 grains with (l x ) : (l y ) : (l z ) = 1:2:2 and Al = 0.2 nm. The condensation energy and 
the spectroscopic gap are displayed also for grains with equidistant single-electron energy 
levels in the interaction band, the spacing between these levels coinciding with the average 
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value (d) for the corresponding ensemble. The ensemble average for the number of electrons 
in the interaction band, (N), varies from 16 at (l x ) = 1.5 nm to 714 at (l x ) = 5.5 nm. 

Fig. 5. Ensemble average {pa,e), where pA,e is given by Eq. QTTJ) , for 1000 superconducting 
grains (A = 0.22) with (l x ) = 1.8 nm, (l y ) = 3.5 nm, (l y ) = 3.7 nm, and Al = 0.25 nm is 
shown as a function of E in the case when the number of electrons in each grain is even 
[panel(a)] and odd [panel(b)]. Inset: function Pa,e(E), averaged over a set of 20 grains with 
(l x ) = 4 nm, (l y ) = 8 nm, (l y ) = 8 nm, and Al = 0.2 nm; the number of electrons in each 
grain is odd. For comparison, the functions (pA,s{E)), calculated for corresponding normal 
grains (A = 0) are also shown. 
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